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ELASTIC MODEL TRANSITIONS: A HYBRID APPROACH 
UTILIZING QUADRATIC INEQUALITY CONSTRAINED LEAST 
SQUARES (LSQI) AND DIRECT SHAPE MAPPING (DSM) 
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A method for transitioning linear time invariant (LTI) models in time varying 
simulation is proposed that utilizes both quadratically constrained least squares 
(LSQI) and Direct Shape Mapping (DSM) algorithms to determine physical dis- 
placements. This approach is applicable to the simulation of the elastic behav- 
ior of launch vehicles and other structures that utilize multiple LTI finite element 
model (FEM) derived mode sets that are propagated throughout time. The time 
invariant nature of the elastic data for discrete segments of the launch vehicle tra- 
jectory presents a problem of how to properly transition between models while 
preserving motion across the transition. In addition, energy may vary between 
flex models when using a truncated mode set. The LSQI-DSM algorithm can ac- 
commodate significant changes in energy between FEM models and carries elastic 
motion across FEM model transitions. Compared with previous approaches, the 
LSQI-DSM algorithm shows improvements ranging from a significant reduction 
to a complete removal of transients across FEM model transitions as well as main- 
taining elastic motion from the prior state. 


INTRODUCTION 

The time invariant nature of elastic data in launch vehicle simulations poses the problem of how 
to transition linear time invariant (LTI) finite element models (FEM) during simulation while pre- 
serving elastic motion across model transitions and ensuring there is no truncation or excitation of 
elastic motion. Launch vehicle controller design must consider the elastic behavior of the system as 
it directly impacts sensor feedback and in turn controller performance. To ensure robust controller 
design, Monte Carlo simulation is performed in both the frequency and time domain to evaluate 
effects of the elastic properties on the system. Simulation of launch vehicle elastic behavior is tradi- 
tionally accomplished using discrete FEM models provided at set intervals throughout the trajectory. 
These models correspond to vehicle mass and other physical properties at given flight conditions 
along the trajectory. Dissimilarities between FEM models requires a transition methodology that 
captures existing elastic motion and does not introduce non-physical effects to the system. 

Previous approaches to transitioning between FEM models relied on linearly interpolating the 
mode eigenvectors; however, there is no mathematical basis for this approach. A method for tran- 
sitioning between multiple FEM models is proposed that utilizes a quadratically constrained least 
squares (LSQI) approach in conjunction with Direct Shape Mapping (DSM) to effect a FEM model 
transition that does not artificially introduce or remove energy from the system. 
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LSQI-DSM is based on the LSQI method, 1 which by way of an energy constraint, uses the eigen- 
vectors of the future state d>fc + i and prior physical state <&/, to estimate the generalized coordinates, 
rj, q, that result in overall system energy at or below the prior state. The two key benefits of the 
LSQI method are that there is no requirement that the models be of similar dimension or that the 
eigenvectors be correlated in any way. 1 

The problem being solved is similar in nature to the topic of FEM model updating discussed in 
the area of structural dynamics 2,3 ’ 2 * 4 ’ 5 ’ 6 ; however, the LSQI-DSM method differs in that it is not 
attempting to update an existing model to experimental data, rather it looks to solve the problem 
of updating to an entirely new model that represents the launch vehicle structure at discrete points 
throughout flight. The LSQI-DSM approach presented is applicable to the simulation of the elastic 
behavior of launch vehicles and other structures that utilize multiple LTI FEM model derived mode 
sets that are propagated throughout time. 

BACKGROUND 

The computational burden of propagating a large number of modes in a launch vehicle simu- 
lation is significant so a subset of modes is typically selected to reduce run time. As a result of 
truncating or using an incomplete modal set, energy between FEM models may not be constant and 
could result in transients immediately following a model transition due to poorly conditioned initial 
states. In an effort to abate these transients, an improved methodology was developed based on 
the quadratic inequality least squares approach, but this new method can accommodate significant 
changes in energy across FEM model transitions. In the LSQI-DSM method, it is proposed that 
physical velocities due to elastic behavior be solved for using the previous LSQI algorithm, but 
solve for displacements using a two-step process that independently addresses the homogeneous 
(un-forced) and steady-state (forced) contributions to the elastic displacement. 

Simulation results of the LSQI elastic model transition method show that while the algorithm 
ensures energy is not introduced across FEM model transitions, E^ + \ < £)., transients “ring” 
the system immediately following a transition, as shown in Figures 2a and 3a. While excitation is 
expected from physical changes to the vehicle configuration during flight such as staging and jettison 
events, substantial excitation due to changes in FEM models at times other than these events are 
non-physical artifacts of the LSQI method. An explanation was sought to determine the root cause 
of these transients and how to abate these effects when transitioning between FEM models. Two 
primary reasons were identified: 

1. Running with a truncated mode set, in particular’ a highly truncated set, invalidates the as- 
sumption that energy is approximately constant across transitions. Figure 1 shows the effect 
of mode count on potential energy across multiple finite element model transitions. 

2. Forcing the physical displacement (qk+ i) and velocity (<jk+ 1 ) of the future state to be the 

final displacement (qQ and velocity (qQ of the prior state meanwhile abruptly changing the 
stiffness and damping to that of the future state results in excitation as the system adjusts to 

instantaneous changes in stiffness and damping. 

FORMULATION OF THE THEORY 

Given a priori knowledge of the stiffness (K), damping (C), and eigenvectors ( *T») from the 
launch vehicle FEM model, the goal is to effect a smooth transition between these models and 
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solve for the physical and generalized state at the point of FEM model transition in a manner that 
preserves existing elastic motion but does not truncate or excite the system in any non-physical 
manner. Prior to presenting the mathematical rationale for the LSQI-DSM method, background on 
the LSQI method is presented as it serves as a partial basis for the LSQI-DSM method. 


Equations of Motion 

The undamped motion is described by: 

Mq + Kq = f (1) 

where M is the mass matrix, K is the stiffness matrix, f is a vector of excitation forces and torques, 
and q is a vector of displacements in physical coordinates. To transform from physical to general- 
ized coordinates the following substitutions arc made to Equation (1): 

q = <f>r) q = $ 7 ) (2) 

to transform the second-order equation in terms of generalized coordinates: 

M<l>77 + = f (3) 

Multiplying Equation (3) by the transpose of the mass-normalized eigenvector matrix ( <E >T ) provides: 

+ $ t K$ 77 = 3> T f (4) 


The eigenvector matrix (4?) is simultaneously selected to satisfy d> T M<l> = I, < h T K < h = diag(w^), 
and the orthonormality constraint. Using this selection of the eigenvector matrix allows the mass 
matrix term 4> T M<I> from Equation (4) to be replaced by the identity matrix (I) and a stiffness 
matrix (K) where the diagonal is comprised of the natural frequencies K) of the elastic modes . 7 


3> T M$ = I 
$ T K3> = K 

where 

/ < 0 ... 0 \ 

~ 0 W n 2 • • • 0 

\o 0 ... c 

A constant viscous damping is typically assumed, leading to a diagonal damping matrix (C). 

/ 2(u ni 0 ... 0 \ 

0 2Cuj n2 ... 0 

C — 

\ 0 0 ... 2 C^n m J 


(5) 

( 6 ) 

(V) 


(8) 


On the right hand side of Equation (4) is the generalized force (S) which represents the effective 
loading associated with all forces and moments on the launch vehicle: 


S = T> T f 


(9) 


Given the generalized coordinate transformation, an identity mass matrix, and constant viscous 
damping. Equation (4) can be rewritten using generalized coordinates: 

77 + C 77 + K 77 = S (10) 
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Energy Constraint 


The least squares problem is bounded by a potential energy constraint on the generalized dis- 
placement (r/) and a kinetic energy constraint on the generalized velocity (r/j. Equation (10) is in 
generalized coordinates and therefore the constraints must also be formulated in generalized coor- 
dinates. Total physical displacement can be written as a sum of all modal contributions, 8 where: <!>, 
is the eigenvector and r]i is the generalized displacement for the i th mode: 

OO 

q = ^2®iVi ( 11 ) 

i= 1 

Similarly, the total potential and kinetic energies are found by summing contributions from all 
modes. Using the transformations from Equation (2), an identity mass matrix (Equation (5)), and 
diagonal stiffness matrix (Equation (7)), the energy equations can be written in physical and gener- 
alized coordinates. 

- OO OO 

C/=-^K^qf = (12) 

i i 

- OO OO 

(12) 

i i 

Equations (12) and (13) represent total kinetic and potential energy assuming an infinitely large 
mode set. For real world applications, the maximum number of modes is finite and dictated by the 
FEM model. As stated earlier, computational requirements necessitate further modal reduction to a 
subset of the complete FEM model. For simulation purposes, the maximum number of modes, m, is 
a finite quantity for which the corresponding “total” potential and kinetic energies are in reality only 
a portion of the total physical quantities. Hereinafter, use of “total kinetic” and “total potential” will 
refer to the total potential ((/') and kinetic (T 1 ) prescribed by the truncated mode set. 


1 m 


Q 2 

i,i K ±i 


= K u*7 i 


(14) 


1 m 1 m 


2 ^ 

i 


(15) 


Equations (14) and (15) can be written as a quadratic form. In addition, since energy is being used 
as a constraint to compare energy both before and after FEM model transitions, the 1/2 multiplier 
can be dropped. Applying these simplifications, the energy Equations (14) and (15) reduce to: 


U' = 77 t Kt7 (16) 

T' = i) J i 7 (17) 


QUADRATICALLY CONSTRAINED LEAST SQUARES (LSQI) METHOD 

The basis functions between FEM models are not equal (4>k / 4?k+iX therefore a method is 
required to calculate the displacement (rjk+i) and velocity (i)k+ i) such that the elastic motion carries 
through the transition and there is no truncation or excitation of the system. Previously the problem 



was formulated as two least squares problems. To bound the solution, constraints on the solution 
limit the resulting magnitude of the kinetic and potential energy. 


minimize||$ fc+ i? 7 fc + i 

Vk+l 

- qfclh 

subject to 

lK + ill 

VI 

(18) 

minimize ||$ fc+ i rj k+1 

Vk+l 

- q/db 

subject to 

ll^+ill 

< wnw 

(19) 


Drawbacks 

The nominal LSQI approach has two primary drawbacks that inhibit its ability to effect smooth 
transitions. The first is that the physically realizable energy should be either equal to or less than the 
energy from the previous state; however, as is shown in the next section, the validity of this assump- 
tion is quite sensitive to the nature of the structure and the number of retained modes. Secondly, it 
does not address the homogeneous and particular solutions to the displacement independently. As a 
result of these issues there is potential for substantial transients across FEM model transitions, as is 
shown in Figures 2a and 3a. 

Non-constant Energy 

The basis behind the quadratically constrained approach was to enforce total system energy con- 
stant across finite element model transitions. Further investigation revealed dramatic changes in 
energy across transitions, especially at low mode counts seen in highly truncated FEM models. 
This is of significant importance because running with a large number of modes adversely impacts 
simulation run-time performance and it is therefore desirable to run with a reduced mode set. To 
better understand the discrepancy in energy between FEM models as well as the number of selected 
modes, the steady-state energy of the system is plotted in Figure 1 for a representative launch vehicle 
to visualize effects of mode count on energy across ten (10) discrete FEM models. 

To estimate steady-state potential, it is assumed the structure is subject to a constant external 
force that has driven the system to an equilibrium state, which is a function of stiffness (K) and 
generalized force (E). Given the potential energy in Equation (16) is a function of generalized 
displacement ( r ] ), the steady-state portion of generalized displacement (r/. s , s ) is found by dropping 
the acceleration (fy) and velocity (ry) terms from Equation (20). 

+ Cr}^ s + Kr] ss = 3>f Kr) ss = E (20) 

Vss = K _1 E (21) 

The steady-state solution ( r ] ss ) represents the particular solution (assuming constant forcing) to 
Equation (20), and since the stiffness (K) is known a priori and if a constant generalized force (E) 
is assumed, then the steady-state potential energy (U' ss ) can be calculated. 

U' ss = vlsK-V ss ( 22 ) 

Previous applications to launch vehicle flex modeling assumed sorting the modes based on modal 
gain placed a majority of the system energy in the first 5% to 10% of modes. However, as shown in 
Figure 1, this assumption does not hold, in particular when using significantly truncated mode sets. 
As mode count increases there is a decrease in A U' between mode sets, as well as a smoothing of 
energy between models. FEM models are unique to their specific application and the method used 
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#Modes #Modes 



Figure 1: Equilibrium potential energy for varying mode counts across multiple FEM models 


for sorting and selecting the mode set directly impacts the elastic response seen in simulation. Given 
the specificity of FEM models to the actual launch vehicle configuration, the analyst must visually 
inspect the energy discrepancy between number of selected modes and manually determine the point 
at which energy convergence is deemed suitable. If constant energy between models is desired, 
a potential method for mode selection could be based on sorting the modes based on equilibrium 
potential and selecting the mode set that ensures energy remains relatively constant between models. 
The drawback of this approach is that the time domain simulation must be capable of supporting a 
variable number of modes throughout the simulation and energy invariance is potentially prioritized 
over other important metrics, such as sensor gain. 

To smooth transitions between FEM models a methodology is needed that addresses the short- 
comings of previous methods by effectively handling abrupt changes in energy: 

■rilKkVk + ril +1 K k+1 rj k+1 (23) 

V J k Vk + vl+^k+i (24) 

as well as considering the homogeneous and steady-state contributions to the total displacement. 

V = Vh + Vss (25) 

LSQI & DIRECT SHAPE MAPPING (DSM) METHOD 

To simplify the formulation, assumptions were made regarding mass matrix normalization (Equa- 
tion (5)), constant external forcing (E/,. = Efc+i), and constant viscous damping (Equation (8)). 
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These simplifications allow for direct solutions to the steady-state displacements through Direct 
Shape Mapping (DSM), which solves for the displacements using the natural frequencies of the 
elastic modes and the external forcing. For the un-forced contributions to displacement, a quadrati- 
cally constrained least squares (LSQI) problem with an energy constraint is solved. 

Steady-State Displacement. Thrust and aerodynamic forces constitute a substantial portion of 
the total launch vehicle system energy and appeal - primarily as potential energy in the system, and 
greatly exceed the kinetic energy in the structure, U' 3> T' . Equation (20) presents the simplifi- 
cations exercised to compute the steady-state portion of generalized displacement (r] ss ), which is 
strictly a function of modal stiffness (K) and generalized force (S), as shown in Equation (21). 
Since the external force is substantial compared to other forces on the launch vehicle, it is important 
to accurately calculate ( rj ss ) to properly capture the forcing contribution. Given the modal frequen- 
cies (u%) are known a priori and the generalized force (S) is assumed constant across transitions 
(Sfc = 3fc + i), it is possible to solve directly for the steady-state portion of the displacement (r/ ss ) 
given the future FEM model. 

Vssk+i = K fc+1 c,k = K fc+1 ($fc + ifk) (26) 

This approach is labeled Direct Shape Matching (DSM) as the generalized displacement following 
the transition (r] ss ) is directly mapped to the equilibrium displacement. 

Homogeneous Displacement. For the homogeneous solution. Equation (10) is set equal to zero: 

Vh + C rj h + K rj h = 0 (27) 

Contrary to the steady-state solution which is a function of the generalized force ( 3 ), the homo- 
geneous solution is a function solely of the damping (C) and stiffness (K). Assuming ( < 1 the 
homogeneous solution can be approximated as a function of only stiffness. To form the constraints 
for the least squares problem, the steady-state displacement from the initial FEM model (r] SSk ) is 
needed to find the homogeneous displacement (r//, ; ) from the same model. 

r] SSk = K~ l S k = K“ 1 (^ fe f k ) (28) 

rth k = Vk ~ Vss k (29) 

The homogeneous displacement is then used to calculate the potential energy due to the homoge- 
neous portion of the solution, which is then used as a constraint when solving for . 

U 'h k = V J hk K k rih k K k+1 = V J hk+1 K k+ ir]h k+1 
The least squares problem can now be solved using the potential energy constraints. 

mininrize||$ fc +i» 7 h fc+1 - ( q * - q SSfc )|| 2 subject to \\U' hk+i || < \\U' hk \\ (30) 

r > h k + 1 

Total Displacement. The total displacement r //. + 1 can now be found as the superposition of the 
steady-state displacement (r] ss ) and the homogeneous displacement (r/h ): 

Vk+1 = Vh k+1 + Vss k+1 (31) 
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Velocity. When solving for the generalized velocity (rjk+i) the external forcing (S) is ignored 
because a constant force only directly impacts displacement ( r /). Since external forcing only directly 
impacts displacement (q) and potential ({/'), the kinetic energy (T 1 ) can be solved for without 
consideration of the generalized force (S). 

T k = V lm U+i = V J k+1 Vk+i 

In a manner identical to the LSQI method, the generalized velocity (r]k+\) is now solved for as a 
least squares problem with a kinetic energy constraint. 

minimize 1 1 fbk+iWk+i — Qfe II2 subject to WU+iW < \\TL\\ (32) 

Vk+l 

ANALYSIS OF LSQI-DSM METHOD 

LSQI-DSM was analyzed using NASA's Space Launch System (SLS) vehicle time domain simu- 
lation tool, MAVERIC (Marshall Aerospace Vehicle Representation in C). To ensure that transients 
at FEM model transitions were reduced, energy was neither added or removed from that dictated 
by each FEM model. The effects of excitation at vehicle transitions on the control system perfor- 
mance were analyzed, and test cases were developed to excite the system over multiple FEM model 
transitions so as to evaluate the effectiveness of the LSQI-DSM method. 

Removal of transients 

In this application, a standalone LSQI approach produces appreciable transients following FEM 
model transitions. Truncation of the mode set invalidates the constant energy assumption between 
models and therefore methods expecting constant energy will not produce smooth transitions. In ad- 
dition, omission of the forced response in the solution to Equation (10) results in poorly conditioned 
initial states at the start of each FEM model, producing severe excitation as shown in Figures 2a 
and 3a, which display the system total energy and sensed acceleration due to elastic motion. The 
data reflects a portion of total flight time and details the first five (5) transitions in a simulation. 




(a) Nominal LSQI (b) LSQI-DSM 


Figure 2: LSQI versus LSQI-DSM and the effect on energy following FEM model transitions 




(a) Nominal LSQI (b) LSQI-DSM 

Figure 3: LSQI versus LSQI-DSM and the effects on sensed acceleration 

In comparison, Figures 2b and 3b highlight the results using the LSQI-DSM method which shows 
improvements ranging from a significant reduction to a complete removal of transients, except for 
cases involving external excitation. 

Verification of energy conservation 

Verification of the LSQI-DSM method was performed by examining elastic motion in the vicinity 
of FEM model transitions. The test case used a sinusoidal force tuned to the first three bending 
mode frequencies and was applied to the system via a force input grid on the vehicle. Intended to 
provide sufficient excitation to oscillate the structure over multiple model transitions, the goal was to 
evaluate LSQI-DSM performance without controller interaction. The test case provides comparison 
of pure DSM, LSQI-DSM, and free decay of the system after being subjected to the sinusoidal 
input. Transitions were disabled to allow estimation of free-decay in the system, while the DSM 
and LSQI-DSM cases were allowed to perform transitions as designed. 

The free-decay case is intended to serve as a benchmark that demonstrates elastic behavior in the 
absence of FEM model model transitions. Figure 4a details the potential and kinetic energies as 
defined by Equations (16) and (17), and highlights the dominance of potential energy in the overall 
system. In the vicinity of the sinusoidal test input, kinetic energy shown in Figure 4b exhibits 
a small amplitude decrease for both DSM and LSQI-DSM; however, LSQI-DSM sustains more 
residual kinetic energy across the transition. For potential energy. Figure 4c details the response for 
all three cases, and Figure 4d provides additional evidence of LSQI-DSM’s ability to sustain energy 
across multiple transitions. For all tests performed, the input force free decay time is consistent. 

Immediately following FEM model transitions, LSQI-DSM sustains kinetic energy across multi- 
ple transitions and shows a decay rate equivalent to that of the free-decay. Preservation of the elastic 
motion and natural decay of the motion is validation that LSQI-DSM is not artificially adding or 
removing kinetic energy. Investigation of the potential energy provides similar conclusions, and 
shows that smooth transition of the potential energy is maintained and is still present two transitions 
following the excitation, meanwhile the DSM method truncates the potential energy. 
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Figure 4: Comparison of various methods and the preservation of elastic motion across transitions 


Impact to vehicle control design 

The SLS control architecture uses an adaptive controller 9 to reduce total controller gain if unde- 
sirable control-structure interactions are detected. With the previous LSQI method, elastic dynamics 
induced at FEM model transitions were sensed by the adaptive controller and interpreted as struc- 
tural instability, leading to a decrease in loop gain. This behavior is undesirable as excitation at 
the model transition is artificial. The LSQI-DSM method reduces/completely removes the artificial 
transient and in turn reduces the controller response to the excitation, as shown in Figures 5a and 5b. 

A drawback of LSQI-DSM when compared to the previous LSQI implementation is discontinuity 
in the physical displacements immediately following a transition. Discontinuity in physical states 
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(a) Nominal LSQI (b) LSQI-DSM 

Figure 5: Positive impact on system total loop gain 


coupled with the vehicle rigid body motion produces a discontinuous kinematic solution, which 
in the presence of inertial navigation, can produce large inertial errors. To solve this problem, 
consideration must be given in the simulation to ensure kinematic continuity across transitions. 

CONCLUSION 

In order to aid in the time domain evaluation, analysis, and verification of launch vehicle control 
system design, a technique to transition between multiple LTI FEM models has been demonstrated. 
The proposed LSQI-DSM method offers appreciable improvements over prior methods by address- 
ing non-constant energy that is experienced when employing a truncated mode set in simulations 
intended to stress specific modes of importance to the control design activities. The LSQI-DSM 
approach accommodates large changes in system energy between FEM models perpetuated by the 
use of a truncated mode set, and shows significant reduction or complete removal of non-physical 
transients, as well as sustaining elastic motion across FEM model transitions. 
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NOTATION 

]jt any variable where subscript indicates the current finite element model 
]fc + i any variable where subscript “k+i” indicates the subsequent finite element model 
] ss steady-state component of the parent variable 
]h homogeneous component of the parent variable 
E total energy 
U potential energy 
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T kinetic energy 

U' potential energy using truncated mode set 
T' kinetic energy using truncated mode set 
q physical displacement 
q physical velocity 
q physical acceleration 
rj generalized displacement 
t] generalized velocity 
i) generalized acceleration 
I identity matrix 
M mass matrix 
K stiffness matrix 
K diagonal stiffness matrix 
C damping matrix 
f excitation forces and torques 
3 generalized force 
$ eigenvectors of the elastic modes 
i mode number 
m maximum number of modes 

natural frequencies of the of elastic modes 
f damping ratio 

REFERENCES 

[1] J. Orr, “Elastic Model Transitions Using Quadratic Inequality Constrained Least Squares,” AIAA Atmo- 
spheric Flight Mechanics Conference, Toronto, CA, 2010. AIAA-2010-7642. 

[2] J. Zapico- Valle, R. Alonso-Camblor, M. Gonzalez-Martinez, and M. Garcla-Dieguez, “A new method 
for finite element model updating in structural dynamics,” Mechanical Systems and Signal Processing, 
2010, pp. 2137-2159. 

[3] J. Mottershead and M. Friswell, “Model Updating in Structural Dynamics: A Survey,” Journal of Sound 
and Vibration 167(2), 1993, pp. 347-375. 

[4] A. Berman and E. Nagy, “Improvements of a Large Analytical Model Using Test Data,” Proceedings of 
the 23rd AIAA/ ASM E/AS C E/AH S Structures, Structural Dynamics and Materials Conference, May 1982. 
Paper 82-0743. 

[5] Y. Yang, Y. Chen, and T. Hsu, “Direct Updating Method for Sttuctural Models Based on Orthogonality 
Constraints,” Mechanics of Advanced Materials and Structures, Taylor & Francis, 2009, pp. 390-401. 

[6] B. Jaishi and W. Ren, “Damage detection by finite element model updating using modal flexibility resid- 
ual,” Journal of Sound and Vibration 290, 2006, pp. 347-375. 

[7] A. Kabe, “Mode Shape Identification and Orthogoanlization,” Proceedings of the 29th 
AIAA/ASME/ASCE/AHS Structures, Structural Dynamics and Materials Conference, May 1988. 
Paper 88-2354. 

[8] D. Hodges and G. Pierce, Introduction to Structural Dynamics and Aeroelasticity. New York, NY: Cam- 
brudge University Press, 2005. 

[9] J. Orr and T. VanZwieten, “Robust, Practical Adaptive ConUol for Launch Vehicles,” August 2012. 
AIAA-2012- 15662. 


12 



Space Launch System 



Elastic Model Transitions: 

^^AJHybrid Approach Utilizing Quadratic Inequality Constrained 
Least Squares (LSQI) and Direct Shape Matching (DSM) 


Robert Jurenko (Leidos Inc. / Jacobs ESSSA) 
Jason Bush (Tri-Vector Services Inc. / Jacobs ESSSA) 
John Ottander (Dynamic Concepts Inc. / Jacobs ESSSA) 


AAS 14-304 


December 16, 2013 


Space Launch System 


Outline 



► Summary 

• Simulation of launch vehicle elastic behavior 

• Previous approaches to transitioning between FEM models 

• Drawbacks and identification of solution path 

► Formulation of the Theory 

• Equations of Motion in generalized coordinates 

• Formulation of energy constraints 

► Nominal LSQI Method 

• LSQI method drawbacks 

► Non-Constant Energy 

► LSQI-DSM Method 

• Steady-state solution 

• Homogeneous solution 

• Total displacement 

► Results 

• Analysis of LSQI-DSM method 

• Removal of transients 

• Verification of energy conservation 

• Impact to vehicle control design 

► Conclusion 

► Acknowledgments 


Elastic Model Transitions: A Hybrid Approach Utilizing Quadratic Inequality Constrained Least Squares (LSQI) and Direct Shape Mapping (DSM) 


2/35 


Space Launch System 


Summary 



The time invariant nature of elastic data in launch vehicle 
simulations poses the problem: 


How to properly transition elastic states between finite element 
models while preserving motion across the transition 

► A method for transitioning linear time invariant (LTI) models in 
time varying simulation is proposed that utilizes: 

• Quadratically constrained least squares (LSQI) and 

• Direct Shape Mapping (DSM) 

► LSQI-DSM applies to the simulation of the elastic behavior of 
launch vehicles and other structures that utilize multiple discrete 
finite element model (FEM) derived mode sets 

► LSQI-DSM algorithm designed to handle truncated mode sets: 

• Can accommodate significant changes in energy across FEM models 

• Ensures flex motion across model transitions 


Elastic Model Transitions: A Hybrid Approach Utilizing Quadratic Inequality Constrained Least Squares (LSQI) and Direct Shape Mapping (DSM) 
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Space Launch System 


Background 



► Launch vehicle controller design considers the elastic 
behavior of the system as it directly impacts sensor 
feedback and controller performance 

► Simulation of launch vehicle elastic behavior is traditionally 
accomplished using multiple FEM models: 

• Models provided at set intervals throughout flight 

• Models correspond to vehicle mass & stiffness 

• Typically provided every 10 to 20 seconds of flight time, for which 
the models are held fixed until conditions reached for next model 

► Dissimilarities between FEM models requires a transition 
method that: 

• Captures existing flex motion of the system 

• Does not introduce non-physical displacement, velocity, or 
acceleration to the system (i.e. ringing) 


Elastic Model Transitions: A Hybrid Approach Utilizing Quadratic Inequality Constrained Least Squares (LSQI) and Direct Shape Mapping (DSM) 
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Background 



► Computational burden of propagating a large number of 
modes in a launch vehicle simulation is significant: 

• Subset of modes is typically selected to reduce run time 

► Using a truncated FEM model can result in: 

• Non-constant energy between models 

• Transients immediately following model transition due to poorly 
conditioned initial conditions 

► The problem to be solved: 

• Given a known set of physical bending displacements and 
velocities at time t 0 , what set of new generalized coordinates 
optimally approximates the physical coordinates using a new 
basis at ti [1] 
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Space Launch System 


Background 



► In search of a physically realistic way to transition between 
models, an approach based on Quadratic Inequality 
Constrained Least Squares (LSQI) [1] was developed to 
effect a smooth transition between models 

• By way of an energy constraint, the LSQI algorithm uses the 
eigenvectors of the future state (<E> fc+ i) and prior physical state 
(<f>fc) to estimate generalized coordinates ( r , /, i)) that result in an 
overall system energy at or below the prior energy level 

• Key benefits of this method include no requirement that the 
models be of similar dimension, or that the eigenvectors be 
correlated in any way 
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l, X (normalized) 


Space Launch System 


Background 



► Simulation results show that while LSQI ensures constant 
energy across a transition, U' k = U k+1 , non-physical 
transients “ring” the system at FEM model transitions 




Time (normalized) Time (normalized) 


Figure 1: Nominal LSQI 


Figure 2: LSQI-DSM 
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Background 



► An explanation was sought to determine the cause of the 
transients at transition points 

► Two primary causes were identified: 

1 . Running with a truncated mode set invalidates the assumption 
that energy is approximately constant across transitions 

•nl^-kVk ^ vl+i^-k+ivk+i (i) 

vlvk vl+iVk +1 (2) 

2. Fixing the physical displacement and velocity across transitions 
while abruptly changing stiffness & damping excites the system 


minimize ^fc+iTjfc+i - 2 s.t. 

Il^+ill <11^11 

(3) 

^fc + l 


minimize — % 2 s.t. 

||2fc+i|| < ||Tfc|| 

(4) 


^fc+l 
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Formulation of the Theory 



► The size and complexity of complete FEM models requires a 
computationally efficient transition methodology 

• The LSQI approach[1] and LSQI-DSM represent efficient 
approaches to updating the FEM model data throughout the 
launch vehicle time-domain simulation 

► Given a priori knowledge of the stiffness (K), damping (C), 
and eigenvectors (<E>) from the launch vehicle finite element 
model, the goal is: 

• To effect a smooth transition between these models 

• Solve for physical/generalized states following model transition 

• Preserves elastic motion across the transition 
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► The undamped motion is described by: 

Mq + Kq = f (5) 

► where M is the mass matrix, K is the stiffness matrix, f is a 
vector of excitation forces and torques, and q is a vector of 
displacements in physical coordinates 

► To transform from physical to generalized coordinates the 
following substitutions are made: 

q = 3?T7 q = 3>r) (6) 

► Equation (5) can now be written in generalized coordinates 

M$rj + K$77 = f (7) 
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Formulation of the Theory 



► Multiplying by the transpose provides: 


$ T M$r7 + $ t K$ 77 = $ T f (8) 


► The eigenvector matrix ($) is selected to satisfy: 

• $ T M$ = I 

• <& T K$ = diag(u;^) 

• and the orthonormality constraint 

► Constant viscous damping is assumed 

0 \ 
o 


K = 


/ “nj 0 

0 WjL 


V 0 


0 \ 

0 


c = 


/ 2Cu n 
0 


V 0 


0 

2 £uj n2 
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► The generalized force represents effective loading of all 
forces & moments on the launch vehicle: 


► Given transformation (6) and prior assumptions, Equation 8 
can be rewritten as: 


► The least squares problem is bounded by: 

• Potential energy constraint on the generalized displacement (77) 

• Kinetic energy constraint on the generalized velocity (77) 

► Recap, total displacement can be written as sum of all 
modal contributions 


= $ T f 


( 9 ) 


77 + C77 + K77 


( 10 ) 


OO 




(ii) 
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Formulation of the Theory 



► Similarly, total potential & kinetic energy are found by 
summing modal contributions 

► Using transformation (6) & earlier assumptions, energy can 
be expressed in generalized coordinates: 

1 OO 1 oo 

C/=-^K M q? = -^K M r 7l 2 (12) 

i i 

1 oo 1 oo 

T = (13) 

i i 

► Equations (12) and (13) represent total kinetic and potential 
energy assuming an infinitely large mode set 

► For simulation purposes, the maximum number of modes, 
to, is a finite quantity 
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Space Launch System 


Formulation of the Theory 



► The corresponding “total” potential and kinetic energies are 
in reality only a portion of the total physical quantities 

► Hereinafter, use of “total kinetic” and “total potential” will 
refer to the total potential ([/') and kinetic (T') prescribed by 
the truncated mode set. 


► Equations (14) and (15) can be written as a quadratic form. 



(14) 


l l 



(15) 


U' = r/ J Kri 
T' = rf r) 


(16) 

(17) 
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Space Launch System 


LSQI Method 



► Basis functions between FEM models not equal 3> k ^ 3>k+i 

► A method is required to calculate the displacement (r) k+ 1 ) 
and velocity (j^+i) given the subsequent FEM model 

► To bound the solution, constraints on the solver limit the 
resulting magnitude of the kinetic and potential energy 


minimize 

\<& k+1 r] k+1 - q fc || 2 

s.t. 

II £4+i II < II *4 II 

(18) 

Uk+l 




minimize 

^fc+i^fc+i — qfc 2 

s.t. 

ll^+ill < ll^ll 

(19) 


Vk+i 
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Energy (normalized) 


Space Launch System 


LSQI Method Drawbacks 



► The original LSQI approach has two key drawbacks that 
generate excitation across FEM model transitions: 

1 . Assumes energy approximately constant across transition 

2. Does not separately address homogeneous & particular solutions 




Figure 3: Nominal LSQI Figure 4: LSQI-DSM 
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Steady- State Potential Energy (n 


Space Launch System 


Non-Constant Energy 



► Investigation of FEM models revealed dramatic changes in 
energy across transitions 



123456789 10 

Finite Element Model # 


Figure 5: Equilibrium potential energy for varying mode counts 
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Non-Constant Energy 



► Previous applications to launch vehicle flex modeling 
assumed sorting the modes based on modal gain 

► FEM models are unique to their specific applications, 
therefore the analyst must visually inspect and manually 
determine the point of energy convergence 

► To provide smooth transitions, a methodology is needed 
that can handle abrupt changes in energy: 

r] T k K k r] k ± ril + 1 K k+ 1 r ) k+1 (20) 

vim ± lil+iVk+i (21) 

► This method must also address both the homogeneous & 
steady-state solutions to the total displacement 

V = Vh + Vss (22) 
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Space Launch System 


Direct Shape Mapping (DSM) 



► Thrust & aerodynamic forces constitute a substantial 
portion of system energy 

► These forces appear in the potential energy and greatly 
outweigh the kinetic energy of the structure, U' > T 

► Given the potential energy in Equation (16) is a function of 
generalized displacement ( 77 ), the steady-state contribution 
of generalized displacement ( r \ ss ) is found by dropping the 
acceleration and velocity terms: 

r^ s + Cp/ + Kr] ss = $f ->• Krj ss = S (23) 
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Space Launch System 


Direct Shape Mapping (DSM) 



► The steadi-state displacement r] ss is strictly a function of 
stiffness (K) and force (E) 

► Given the modal frequencies “n are known a priori and the 
generalized force E is assumed constant across transitions: 

Efc = 3 fc+1 (24) 

steady-state displacement can be solved directly: 

Vss k+ 1 = K fe+1 £ifc (25) 

► This this approach is labeled Direct Shape Matching (DSM) 
as the generalized displacements are directly mapped to the 
equilibrium displacement governed by ® k +i, w nfc+1 , and E fe 
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Space Launch System 


LSQI-DSM 



► Homogeneous solution -> Equation (10) set equal to zero: 

• Function solely of damping (C) and stiffness (K). 

• Assuming ( < 1 the homogeneous solution can be approximated 
as a function of only stiffness 

Vh + Crj h + K r] h = 0 (26) 

► To form the constraints for LSQI, the steady-state 
displacement from the initial FEM model (r] ss k ) is needed to 
find the homogeneous portion of the displacement (rj hk ) 

rj SSk =K- 1 H fe = K- 1 ($ fc f,) (27) 

Vh k =rj k - r] SSk (28) 
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Space Launch System 


LSQI-DSM 



► The homogeneous displacement is used to calculate the 
homogeneous portion of potential energy 


U L =r il k ^kVh k 

( 29 ) 

^ h k +i Vh k+1 '^-k+l'nh k+1 

( 30 ) 


► The least squares problem is now solved subject to the 
energy constraints: 

minimize||# fc+ i»7 hfc+1 -(q fc -q SSfe )|| 2 s.t. \\U' h || < \\U' h J (31) 

r,h k + 1 

► Total displacement rj is the superposition of steady-state 
and homogeneous solutions: 

Vk+1 = T]h k+1 + r] SSk+1 ( 32 ) 
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Space Launch System 


LSQI-DSM 



► When solving for generalized velocity (rj k+ 1) the external 
forcing is ignored 

• Constant forcing only directly impacts displacement (rj) 

• Kinetic energy solved without consideration of generalized force 

T k = Vim T' k+1 = vl+iVk+i (33) 

► Similar to nominal LSQI, generalized velocity is solved as a 
least squares problem with a kinetic energy constraint 

minimizej|$ fc+1 7) fe+1 -q fc || 2 s.t. ||l£ +1 || < \\r k \\ (34) 

m + 1 
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SIS Analysis of LSQI-DSM Method 

Space Launch System 


► LSQI-DSM was evaluated using MAVERIC (Marshall 
Aerospace Vehicle Representation in C) 

► Two problems drove excitation in LSQI approach: 

• Truncated mode set invalidates constant energy assumption 

• Omission of forced response in the displacement solution resulted 
in poorly conditioned initial states at the start of each FEM model 

► LSQI-DSM method addressed these two issues by: 

• Using DSM to initialize the forced response 

• Using the LSQI approach to capture the un-forced response 

• Total solution is a combination of forced & un-forced solutions 
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Energy (normalized) 


Space Launch System 


Removal of Transients 



► LSQI-DSM shows improvements ranging from a significant 
reduction to a complete removal of transients, except for 
cases involving external excitation 

► In case of energy, LSQI-DSM transitions without excitation 




Figure 6: LSOI Figure 7: LSQI-DSM 
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l, X (normalized) 


Space Launch System 


Removal of Transients 



► In case of sensed acceleration, LSQI-DSM significantly 
reduces excitation, which has: 

• Positive impact on sensor output 

• Leads to reduced error in navigation solution 



Time (normalized) 



Figure 8: LSQI 


Figure 9: LSQI-DSM 
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SIS Verification of Energy Conservation 

Space Launch System 


► Verification of LSQI-DSM approach was performed by 
examining elastic motion around FEM model transitions 

► Test case: 

• Sinusoidal force tuned to the first three bending mode frequencies 

• Force applied via an input grid on the vehicle and intended to 
provide sufficient excitation over multiple FEM model transitions 

• Goal of test case was to evaluate LSQI-DSM performance without 
controller interaction 

► The test case provides a comparison of pure DSM, 
LSQI-DSM, and free decay of the system 

► Free-decay case used as benchmark to demonstrate elastic 
behavior in the absence of FEM model transitions 
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Energy (Potential and Kinetic) (normalized) 


Verification of Energy Conservation 



Space Launch System 


► Immediately following FEM model transitions: 


• Kinetic energy exhibits small amplitude decrease for both methods 

• LSQI-DSM sustains more residual energy across the transition 




Figure 10: LSQI 


Figure 11: LSQI-DSM 
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Potential Energy (normalized) 


Verification of Energy Conservation 



Space Launch System 


• LSQI-DSM sustains energy across multiple transitions 

• Shows a decay rate equivalent to that of the tree-decay 

• These facts provide validation that LSQI-DSM is not artificially 
adding or removing energy 


||||||^ — ~ 


Free Decay 

DSM 

- - LSQI-DSM 



Time (normalized) 

Figure 12: LSQ1 
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SIS Impact to Vehicle Control Design 44 

Space Launch System 


► SLS uses an adaptive controller to modify gains[2]: 

• Adaptive controller interprets transients as structural instability 

• LSQI-DSM reduces/completely removes the artificial transients, 
effectively reducing the controller response to the excitation 




Figure 14: LSOI 


Figure 15: LSQI-DSM 
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Space Launch System 


Drawbacks to LSQI-DSM 



► Compared to previous LSQI implementation, LSQI-DSM 
produces discontinuity in the physical displacements 
immediately following a transition: 

1 . Results in discontinuous kinematic solution 

2. In presence of inertial navigation, this discontinuity can lead to 
large inertial navigation errors 

► To solve this problem, consideration in the simulation must 
be given to ensure kinematic continuity across transitions 
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Space Launch System 


Conclusion 



► The LSQI-DSM method offers appreciable improvement over 
prior methods by: 

• Accommodating large changes in system energy between FEM 
models perpetuated by the use of a truncated mode set 

• Independently addressing the unforced and forced components of 
the displacement 

► LSQI-DSM shows significant reduction or complete removal 
of non-physical transients following FEM model transitions: 

• Positive impact to controller response and navigation solution 

• Sustains elastic motion across multiple FEM model transitions 
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Backup 



► The following slides are backup 
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